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LONG-TERM  GOALS 

Terrestrial  runoff  and  river  input  dominates  urban  pollutant  loading  rates  degrading  nearshore  and 
surfzone  water  quality  ( e.g Boehm  el  al.,  2002).  Surfzone  mixing  processes  disperse  and  di¬ 
lute  such  (and  other  types  of)  pollution.  On  smaller  length-scales  (smaller  than  the  water  depth), 
breaking-waves  and  bed-generated  turbulence  mix  tracer.  However,  field  surfzone  observations  of 
turbulence  previously  have  been  extremely  scarce,  and  much  about  surfzone  small-scale  turbulence 
is  not  known.  On  larger  scales  (10-100  m),  horizontal  dispersion  is  driven  by  surfzone  eddies  and 
meanders  associated  with  shear  waves  ( Oltman-Shay  et  al.,  1989)  or  finite  breaking  crest  length 
(. Peregrine ,  1998).  Understanding  the  small  and  large  length-scale  mixing  processes  important  to 
predicting  the  fate  (transport,  dispersal,  and  dilution)  of  surfzone  tracers  whether  pollution,  bacte¬ 
ria,  larvae,  or  nutrients. 

OBJECTIVES 

The  scientific  objective  is  to  improve  understanding  and  modeling  of  dispersion  of  tracers  (pol¬ 
lution,  fecal  indicator  bacteria,  fine  sediments)  within  the  nearshore  (a  few  100  m  of  the  shoreline) 
and  the  surfzone.  In  this  report,  the  focus  is  on  three  research  components  built  upon  observations 
from  the  HB06  experiment  (Pis:  Feddersen  and  Guza).  The  first  is  stochastic  modeling  of  surfzone 
drifter  dispersion  from  the  HB06  experiment  ( Spydell  and  Feddersen ,  201 1).  Second,  studying  the 
small-scale  turbulence  in  the  surfzone  due  to  breaking  waves  and  bottom  boundary  layer  processes 
(. Feddersen ,  2011).  Third,  is  modeling  nearshore  nutrient  fluxes  and  the  resulting  phytoplankton 
and  comparing  it  to  observations  ( Omand  et  al.,  2011b).  In  addition,  IB09  experiment  (performed 
in  collaboration  with  R.  T.  Guza)  analysis  is  ongoing  and  is  not  described  here. 

APPROACH 

HB06  Dye  Dispersion  Modeling 

A  Boussinesq  wave-current  model  funwaveC  has  been  coupled  with  a  tracer  evolution  model. 
The  model  reproduces  the  observed  wave  and  current  conditions  on  the  5  days  of  HB06  dye  re¬ 
leases.  Dye  tracer  is  released  in  the  model,  and  the  model  dye  tracer  transport  and  dispersion  is 
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analyzed  analogous  to  the  observations  ( e.g .,  Clark  et  al.,  2010a).  The  results  are  presented  in 
Feddersen  et  al.  (201 1)  and  Clark  et  al.  (201 1),  and  are  not  discussed  further  here. 

Stochastic  Particle  Simulation  for  Surf  zone  Dispersion 

Drifter- derived  diffusivities  are  time-dependent.  In  an  unbounded  domain,  diffusivities  are 
monotonic  with  a  linear  ballistic  regime  for  times  -C  the  Lagrangian  time-scale  r  and  become 
constant  at  times  3>  r.  For  five  HB06  experiment  drifter  release-days  on  an  alongshore  uniform 
beach,  the  cross-shore  diffusivity  Kx  was  estimated  for  times  <  1000  s  using  an  unbiased  estimator 
(, Spydell  et  al.,  2009).  The  estimated  Kx  had  a  ballistic  regime  for  t  <  50  s,  a  maximum  around 
100  s,  with  a  slow  decay  for  longer  times.  A  potential  reasons  why  a  monotonic  Kx  may  not 
be  expected  in  the  surfzone  is  the  presence  of  a  boundary  of  the  beach  which  prevents  unbounded 
diffusion.  Here  HB06  particle  trajectories  are  stochastically  simulated  with  the  Langevin  equations 
with  a  shoreline  boundary  to  explain  the  observed  features  of  the  cross-  and  alongshore  diffusivities 
The  approach  here  was  developed  in  collaboration  with  Project  Scientist  Matt  Spydell. 

Small-scale  Surfzone  Turbulence 

The  vertical  structure  of  turbulence  in  the  surfzone  is  of  interest.  Both  breaking  waves  and 
near-sea-bed  shear  are  possible  sources  of  turbulence.  Here  a  key  turbulence  statistic,  the  turbulent 
dissipation  rate  e  is  estimated  from  Acoustic  Doppler  Velocimeters  observations  following  Fedder¬ 
sen  et  al.  (2007)  and  Feddersen  (2010).  The  effect  of  both  wave -breaking  and  bottom  boundary 
layer  processes  (BBL)  upon  e  are  examined. 

Episodic  Nutrient  Fluxes  and  Phytoplankton  Modeling 

In  Southern  California,  intense  phytoplankton  blooms  localized  in  the  nearshore  (<  20  m 
depth)  appear  intermittently,  particularly  during  summer  and  fall.  The  underlying  drivers  of  these 
blooms  are  poorly  understood.  Three  distinct  phytoplankton  blooms  lasting  4-9  days  were  ob¬ 
served  in  approximately  15  m  water  depth  near  Huntington  Beach  CA  between  June  -  October  of 
2006  during  the  HB06  experiment.  Vertical  nutrient  (nitrate)  fluxes  are  parameterized  and  used  to 
drive  a  nitrate -phytoplankton  model.  This  work  is  part  of  former  graduate  student  Melissa  Omand’s 
thesis. 

WORK  COMPLETED 

•  Clark  et  al.  (2010a)  of  the  HB06  dye  dispersion  studies  has  been  published  in  JGR  Oceans. 

•  A  manuscript  ( Feddersen ,  2010)  has  been  published  in  J.  Atmospheric  and  Oceanic  Tech.. 
This  manuscript  deals  with  the  methods  of  analyzing  Acoustic  Doppler  Velocimeter  data  for 
estimating  the  turbulent  dissipation  rate  in  the  surfzone  and  air-sea  boundary  regions 

•  Omand  et  al.  (201  la)  has  been  published  in  Limnology  and  Oceanography,  reporting  on  the 
evolution  and  dynamics  of  a  nearshore  red  tide  observed  during  HB06. 

•  Feddersen  (2011),  on  surfzone  turbulence  dynamics  from  HB06  observations,  is  press  to 
JPO. 
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•  Spydell  and  Feddersen  (2011),  on  a  theoretical  analysis  of  finite-Lagrangian  time-scale  on 
shear  dispersion  is  in  press  to  JFM. 

•  A  two-part  paper  series  ( Clark  et  al.,  2011;  Feddersen  et  al.,  201 1)  on  Boussinesq  modeling 
of  surfzone  tracer  dispersion  is  in  press  to  JGR. 

•  A  manuscript  ( Spydell  and  Feddersen ,  2011)  on  stochastic  simulation  of  surfzone  particle 
dispersion  is  in  preparation  for  JGR. 

•  Analysis  of  the  IB09  (Imperial  Beach  CA  in  Sept-Oct  2009)  experiment  (Pis:  Guza  and 
Feddersen)  is  in  full  swing. 


RESULTS 

Background:  HB06  Experiment 

Observations  were  collected  from  15  September  to  17  October  2006  (800  hours)  at  Hunting- 
ton  Beach  CA,  a  site  with  chronic  water  quality  problems.  A  cross-shore  transect  of  co-located 
pressure  sensors  and  acoustic  Doppler  Velocimeters  was  deployed  spanning  160  m  out  to  4  m 
mean  water  depth.  The  tide  range  was  nominally  ±1  m.  The  data  was  sampled  at  8  Hz.  The 
ADVs  sampled  between  0.5-1. 0  m  above  the  bed.  The  cross-  and  alongshore  coordinate  are  x 
and  y,  respectively.  The  mean  water  depth  is  given  by  h.  At  each  of  the  frames,  hourly  estimates 
of  significant  wave  height  Hs ig,  mean  alongshore  current  v,  and  turbulent  dissipation  rate  e  were 
estimated. 


Stochastic  Simulation  of  Surfzone  Drifter  Dispersion 

An  unusual  feature  of  the  HB06  cross-shore  drifter  diffusivities  is  their  non-monotonic  nature 
(they  decrease  with  time).  This  behavior  is  explored  by  stochastic  particle  simulations  governed 
by  the  Langevin  equations.  Please  see  Spydell  and  Feddersen  (2011)  for  full  details. 

The  results  of  the  model  simulations  are  shown  in  Figure  1.  At  short  times  (t  <  50  s  and 
t / tx  <  0.7)  the  modeled  Kfu)  and  analytic  Kx‘ ’  diffusivities  reproduce  well  the  observed  Kf  '1 
(compare  the  dashed,  solid,  and  dash-dot  curves  in  Figure  1,  left  column)  as  cross-dispersion  is  in 
a  ballistic  regime  (Kx  =  alt).  This  also  indicates  that  the  bulk  cross-shore  velocity  variance  ax 
used  in  the  modeled  and  analytic  solutions  is  accurate.  With  the  exception  of  10/14,  the  modeled 
/vfni)  reproduces  the  observed  K(f>  at  longer  times  (t  >  400  s  or  t/rx  >  3,  left  and  right  columns 
of  Figure  1,  respectively).  Similarly,  at  longer  times,  the  analytic  also  reproduces  Kf  '1  on 
10/02  and  10/03,  but  somewhat  underpredicts  K^>  on  09/17  and  10/15.  As  discussed  in  Spydell 
et  al.  (2009),  the  poor  model  prediction  on  10/14  at  longer  times,  drifter  trajectories  converged  in 
the  inner-mid  surfzone  suggesting  bathymetric  control  and  resulting  in  the  rapid  longer  time  decay 
in  K(°\  For  more  information  please  see  Spydell  and  Feddersen  (2011). 

Small-scale  Surfzone  Turbulence:  Dissipation  and  Wave  Energy  Flux  Relationship 
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Figure  1 :  (left)  The  cross-shore  diffusivity  Kx  versus  time  t  and  (right)  the  non-dimensionalized 
diffusivity  Kx/kx  versus  t/rx  for  the  5  dye-release  days  (rows).  In  each  panel  the  observed 
(dashed)  and  Langevin  equation  modeled  (solid)  and  the  analytic  (dash-dot)  are  indicated  by  the 
legend.  In  all  panels  the  shaded  gray  region  represents  the  observed  K^>  sampling  error.  In  panels 
(a2)-(e2),  the  best-fit  tx  and  and  kx  =  a^Tx  are  used,  and  the  analytic  expression  for  the  long-time 
Kx/kx  =  (1  —  2 / 7r)  ~  0.36  is  shown  in  the  thin  dashed  line. 
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During  the  HB06  experiment,  the  turbulent  dissipation  rate  e  typically  was  significantly  larger 
(by  a  factor  of  10)  within  the  surfzone  relative  to  seaward  of  the  surfzone,  suggesting  the  impor¬ 
tance  of  wave  breaking  to  surfzone  turbulence.  Waves  approaching  the  surfzone  have  an  associated 
onshore  wave  energy  flux  F  that  is  conserved  until  wave  breaking  begins.  As  F  =  0  at  the  shore¬ 
line,  the  incoming  wave  energy  must  be  converted  into  other  forms  of  energy  within  the  surfzone. 
In  the  simplest  steady-state  energy  balance,  the  incident  wave  energy  flux  is  balanced  by  the  depth- 
integrated  turbulent  dissipation  over  the  entire  surfzone.  If  dissipation  were  depth-uniform,  then 
the  simple  cross-shore  integrated  energy  balance  between  the  incident  wave  energy  flux  and  surf¬ 
zone  dissipation  becomes 

rx  b 

/  hedx  =  FXF7.  (1) 

J  0 

The  dissipation  rate  e  varies  over  the  vertical  (e.g.,  George  et  al.,  1994),  and  surfzone  laboratory 
experiments  indicate  the  majority  of  dissipation  occurs  above  trough  level  as  found  in  laboratory 
studies..  Thus  the  assumption  that  /  e(z)dz  =  he  is  not  appropriate.  However  the  observed  he 
likely  is  proportional  to  the  depth-integrated  dissipation  (e.g.,  Trowbridge  and  Elgar,  2001),  par¬ 
ticularly  as  e  co-varies  across  the  array  (Fig.  2).  Therefore  the  balance  /  he  dx  =  cFx F7  (similar  to 
Eq.  8)  is  examined  where  c  is  a  fit  constant  of  proportionality. 

Due  to  data  gaps,  J(]Ch  he  dx  is  calculated  in  two  manners:  The  first  only  estimates  the  integral 
when  all  surfzone  frames  have  good  e  estimates  resulting  in  IV  =  143  data  points.  The  second 
requires  at  least  2  (for  Xb  at  F3  or  F4)  or  3  (for  Xb  at  F5  or  F6)  good  surfzone  e  to  calculate  the 
integral  resulting  in  N  =  430  data  points. 

The  integrated  surfzone  dissipation  J((h  he  dx  using  either  estimator  is  linearly  related  to  the 
incoming  wave  energy  flux  FXf7  (Fig.  2),  demonstrating  the  link  between  incoming  wave  en¬ 
ergy  and  viscous  dissipation  to  heat,  but  /0Xb  he  dx  is  two  orders  of  magnitude  smaller  than  FXf7 
(Fig.  2).  With  the  first  /0Xb  he  dx  estimate  (N  =  143),  the  relationship  between  FXf7  and  he  dx 
has  moderately  high  squared  correlation  r2  =  0.61  with  least-squares  best-fit  slope  of  c  =  0.01 
(Fig.  2a),  indicating  that  only  1%  of  the  depth-normalized  wave  energy  is  observed.  Using  the 
second  Jgb  he  dx  estimator,  with  3x  the  number  of  good  data  points  (N  =  430),  the  relationship  is 
similar,  but  noisier,  with  squared  correlation  r2  =  0.35  and  slope  of  c  =  0.008  (Fig.  2b).  For  more 
details  please  see  Feddersen  (201 1). 

Episodic  Nutrient  Fluxes  and  Phytoplankton  Modeling 

Here,  we  estimate  the  vertical  advective  and  turbulent  nitrate  fluxes  during  brief  (<  1  week) 
pulses  of  N03  to  the  nearshore  euphoric  zone  in  18  m  water  depth,  between  mid- June  and  mid- 
October  2006,  at  Huntington  Beach,  CA.  A  nearby  Clil  a  timeseries  chronicles  the  rapid  growth 
and  decline  of  three  distinct  phytoplankton  blooms,  one  of  which  was  a  red  tide  of  the  dinoflagel- 
lates  Lingulodinium  polyedrum.  We  use  temperature  as  a  proxy  for  N03,  and  inferred  the  episodic 
turbulent,  and  advective  fluxes  from  water  column  measurements  of  currents  and  temperature. 
Each  N03  flux  event  preceded  a  bloom  event,  indicating  that  blooms  may  be  a  response  to  these 
fluxes.  The  correlation  between  the  N03  flux  and  the  observed  C 1 1 1  a  was  maximum  (r2  =  0.40) 
with  an  8  day  lag.  A  simple  local  Nitrate-Phytoplankton  (NP)  model  using  a  linear  uptake  func¬ 
tion  and  driven  with  the  NO3  flux  captured  the  timing,  magnitude,  and  duration  of  the  three  Chic/ 
blooms  (skill=  0.61)  using  optimal  net  growth  rate  parameters  that  were  within  the  expected  range 
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(Figure  3).  The  success  of  a  very  simple  2  parameter  NP  model  in  reproducing  the  fundamental 
features  of  all  three  blooms  are  included)  highlights  the  strong  connection  between  the  vertical 
nitrate  flux  and  the  lagged  Chi  a  response,  and  may  assist  the  design  of  future  nearshore  programs 
identify  the  critical  physical  parameters  and  timescales  to  gain  a  potentially  predictive  insight  into 
bloom  dynamics  in  Southern  California.  For  more  details  please  see  Omand  et  al.  (201  lb). 

IMPACT/APPLICATIONS 

Potential  impacts  include  improving  surfzone  and  nearshore  mixing  parameterizations  based 
upon  bulk  factors  such  as  wave  height,  wave  period,  bathymetry,  and  currents. 

RELATED  PROJECTS 

The  Tidal-Inlets/River-Mouths  DRI  project  is  building  upon  the  dye  tracer  and  drifter  results 
here,  but  with  expanded  geographical  scope  to  include  a  tidal  inlet. 
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Figure  2:  Cross-shore  integrated  energy  balance:  /0Cb  he  da;  versus  incident  F7  wave  energy  flux 
FXf7  for  fgb  he  da;  calculated  when  (a)  all  surfzone  e  are  good  (N  =  143)  and  (b)  at  least  2  ( xb  at  F3 
or  F4)  or  3  (if  Xb  is  at  F5  or  F6)  surfzone  e  are  good  ( N  =  430).  The  black  dashed  line  represents 
the  least-squares  best-fit  constrained  to  go  through  the  origin  with  the  best-fit  slope  c  and  squared 
correlation  r2  of  (a)  c  =  0.01  and  r2  =  0.61  and  (b)  c  =  0.008  and  r2  =  0.35,  respectively. 
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Figure  3:  Timeseries  of  24  h  averaged  observed  Chla(obs)  (/igL-1,  light  green  line)  and  modeled 
Chi  a(m)  (/ig  L-1,  dark  green  line)  Chlorophyll-a.  The  gray  shaded  regions  indicate  the  time  periods 
of  the  3  bloom  events. 
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